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During this period a paper 1 was presented at the AIAA 10th Applied Aerodynamics 
Meeting. Later, the paper was submitted for publication in the AIAA Journal and was accepted. 
In addition, an abstract 11 was submitted to the 24th AIAA Fluid Dynamics Conference. Copy of 
the abstract is enclosed. 
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Introduction 

It is well known that current k—t models cannot predict the flow over a fiat plate 
and its wake 1-2 . In an effort to address this issue and other issues associated 
with turbulence closure, a new approach for turbulence modeling is proposed 
which exploits similarities in the flow held. Thus, if we consider the flow over 
a flat plate and its wake, then in addition to taking advantage of the log-law 
region, we can exploit the fact that the flow becomes self-similar in the lar 
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wake. This latter behavior makes it possible to cast the governing equations 
as a set of total differential equations. Solutions of this set and comparison 
with measured shear stress and velocity profiles yields the desired set of model 
constants. Such a set is, in general, different from other sets of model constants. 
The rational for such an approach is that if we can correctly model the flow over 
a flat plate and its far wake, then we can have a better chance of predicting the 
behavior in between. It is to be noted that the approach does not appeal, in any 
way, to the decay of homogeneous turbulence. This is because the asymptotic 
behavior of the flow under consideration is not representative of the decay of 
homogeneous turbulence. 


Approach 

The approach will be illustrated by a k model. A Similar approach can be 
used for a k — t model. For a plane wake, the defect velocity can be represented 
as 3 

U e - u = U s f (y) = ( r sf(rj) ( L ) 

where i' e is the free stream velocity, u is the velocity, and / is a length scale 
representative of the width of the wake. For the wake 

U s ~ ~ /2 ( 2 ) 


The k and u> equations can be written as 
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The mean velocity, u is governed by 
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In the above equation, v is the kinematic viscosity. (T m . ( T 1 . C ^ are model 
constants, and the subscript t designates a turbulent quantity. Using Equation 
(1) and setting k and u> as 

k = i' s 2 h(i /), v t = C\- (6) 

the equations in the far wake reduce to a set of second order total differential 
equation in ?/. Two dimensionless parameters appear in these equations u/L'J 
and UJ/UeX. The first parameter was chosen over the range 10~ 4 — 10 -ti The 
second parameter can be expressed as 

UJ _ U, l 0 
T~c ~ T~ s 0x 

= ‘ 2 )jr§ ( 7 ) 

where 9 is the momentum thickness. As was indicated in Ret. [1], the general 
concensus is that v t j U e 9 has a value of .032 when v t is assumed constant. In 
this work v t is variable. Because of this, the quantity v to jU,0, where i/ to is 
the value of v t at the line of symmetry, is assumed to be .032. The boundary 
conditions are 

/'(*/) = h'(v) = </'('/) = o at >1 = 0 (8) 

f(v) — " />('/) -*• g{r]) 0 as t I -* oc (9) 

Results and Discussion 

The above system has four model constants. Because we are interested in model 
constants that are suited for Hows over the body and its wake, the results of the 
log-law region are used to relate some of the model constants, l he resulting 
equations can be written as 

■ ) 

( 10 ) 

where k is the Karman constant and (.\ t = .00. Using Equation! 10). the re- 
maining constants were obtained by finding the "best" set of constants that 
reproduced the shear stress distribution, i.e. 



and the velocity defect distribution, Ur])- Because we are dealing with a two- 
point boundary value problem, the nature of the solutions obtained are shown 
in Figure (1). The selection of the model constants reduces to selecting the 
solution that has the shape labeled "correct solution" and gives a best fit of 
shear and defect velocity distributions. Since e , 0 / f V 0 is assumed to be .032. 
the values of g(0) and h(0) are related by the relation 
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— = C — — not 

U e 0 *l\6g(0) 

The quantity HJ s /U e 0 is obtained from Ref. [1] as [( 4 7r ) ln{2)} 2 

The approach outlined above is used to calculate the model constant in a k— 
uj model. Figure (2) shows the distribution of k . ( /? ( 77 ) ) and uu (g(7/)). The shear 
stress distribution and velocity profiles are compared in figure (3) and (4) with 
the asymptotic solution, where v t is assumed constant throughout, and with 
the experiments of Ramaprian et al. 4 and Pot 0 . A summary of the constants 
obtained is given in Table 1. The table contains Wilcox's constants 0 along 
with constants developed by Anderson'. The results obtained by Anderson 
employed numerical optimization of computations starting on the flat plate 
and extending to | = 2.5xl0 3 downstream of the plate trailing edge, and by 
comparison with asymptotic solution and available experiment. 


Proposed Work 

I he next step in the research is to repeat the procedure for a k — t model. I he 
resulting sets of constants will then be used to calculate the flow held over the 
plate and its wake. Results will be compared with available experiment. 
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Table 1 . Comparison of Various Model Constants 


Constant 

Anderson 

Wilcox 

Current 

Cv! 

0.80 

1.5667 

0.49 1 


1.1733 

1.8333 

0.8633 

Cji 

0.09 

0.09 

0.09 

G k 

1.17647 

2.00 

1.17647 

&UJ 

1.42857 

2.00 

1.42857 
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Figure 1 . Method of Determining Paremeters 
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Figure 2. T.K.E. h (rj) and Frequency Scale g(//) Proiiles 
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Figure 3. Shear Stress Profile. (i(//) 


9 


i oor 

0.75 

m 

0.50 

0.25 

0.00 

0 . 



X. 0 

Hu Ramaprian 99 

+ Ramaprian 23.1 

- Ramaprian 59.3 

X Ramaprian "9.1 

^ : Pot 960 960 

Current 

“ Asymptotic 




0 5 


n=y/i 


1.0 


15 


Figure 4. Velocity Defect Profile. f(//) 
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